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Abstract 

B — T phase diagram with a multicritical point for interlayer Josephson vortices is mapped out based on 
Monte Carlo simulations. For high magnetic fields we find a novel Kosterlitz-Thouless (KT) intermediate phase 
characterized by in-plane two-dimensional (2D), quasi long-range orders (QLROs) of vortex alignment and 
superconductivity. Both the transition to high-temerature normal state and the evolution to low-temperature 
phase of 3D LRO are continuous. Decoupling of the 3D system into the 2D state is triggered by hops of segments 
of Josephson flux lines across superconducting layers activated by thermal fluctuations. For low magnetic fields, 
a single first-order melting transition is observed. 



PACS numbers: 74.60. Ge, 74.20.De, 74.25Bt, 74.25.Dw 



The discovery of high-T c superconductivity in cuprates has triggered extensive researches into the finite- 
magnetic-field superconductivity transition in type 11 superconductors [1]. It is now well established that the 
transition is first order, accompanied by the freezing of the flux-line liquid into lattice. This notion is particularly 
important because the transition was considered second order for a long time since Abrikosov [2] . 

The high-T c superconductors share a layered structure in which the superconductivity is widely believed to 
occur mainly in the CuC>2 layers intervened by layers of charge reservoir. This profound layer structure causes 
significant differences in both equilibrium and transport properties under magnetic fields in different directions. 
The first-order normal to superconductivity transition is mainly observed under magnetic fields perpendicular to 
the layers, at which the two-dimensional (2D) translation symmetry enjoyed by pancake vortices is broken. In a 
sharp contrast, a parallel magnetic field penetrates the system through the reservoir layers, in terms of Josephson 
vortices. The relevant c-axis translation symmetry is broken a priori, which may raise new phase and new melting 
process. In fact, a peculiar transport phenomenon has been found in Bi 2 Sr 2 CaCu20 8+ j / by lye et al. [3] that the 
resistivity does not depend on the angle between the magnetic field and current when they are both parallel to the 
CuC>2 layer [3] and by Ando et al. that the system shows power-law, non-Ohmic dissipations [4]. Chakravarty et al. 
pointed out that such a dissipation cannot occur in a lattice phase [5] . Blatter et al. [6] proposed a novel Kosterlitz- 
Thouless (KT) [7, 8] scenario at high magnetic fields, characterizing the behavior by a smectic state with vanishing 
interlayer shear modulus (see [9] for a possible KT phase at intermediate magnetic fields). Motivated by the 
experiment by Kwok et al. suggestive of continuous melting transition [10], Balents and Nelson proposed a smectic 
phase in between the lattice and liquid states and thus a two-step, continuous melting scenario at intermediate 
magnetic fields [II] (see [12] for a proposal of supersolid). Decoupling between superconducting layers was actually 
proposed even before the discovery of high-T c superconductors [13, 14]. Starting from a 2D elastic Hamiltonian 
and resorting to renormalization group (RG), however, it was argued by Mikheev and Kolomeisky [15] that only 
a 3D long-range crystalline order is possible besides liquid state (see also [16]). The discrepancy among different 
approaches has not been resolved to provide a unified picture, though they provide new physical insights, 

By means of computer simulations we find in the present work that there is a multicritical point in the B — 
T phase diagram of interlayer Josephson vortices: Below the critical field, a single first-order transition upon 
temperature sweeping is observed; above it, there exists an intermediate KT phase characterized by in-plane 2D, 
quasi long-range orders (QLROs) of vortex alignment and superconductivity in between the normal phase and 3D 
lattice phase, accompanied by two continuous melting transitions. 

The Hamiltonian is for the phases of the superconductivity order parameter on the simple cubic lattice [17, 18] 

H = -J c °s(^-<ft)-^ cos(ipi - (fj - ^ / A c dr c ), (1) 

<i,j> ||ic,yaxis <2,j>||caxis 

where J = (/>g<i/167r 3 A^ h , the y direction along the external magnetic field and x _L c _L y, A — (0,0, —xB), and 
7 = A c /A a fc. The system is of size L x x L y x L c = 384d x 200c? x 20d under periodic boundary conditions, with 
d the separation between Cu02 layers. Vortices are identified by counting gauge invariant phase differences [17]. 
It was found [18] that the ((-function peak in the specific heat associated with the first-order melting of Josephson 
vortex lattice is suppressed when the magnetic field and/or anisotropy parameter 7 exceed the critical values given 
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by fj = 1/2a/3 (/ = Bd 2 /4>o), suggestive of a multicritical point and continuous melting(s) beyond it. In the 
present study, we fix / = 1/32 and take two typical anisotropy constants: 7 = 8, realized in YBa2Cu307_«5 and 
slightly below the critical value for / = 1/32, and 7 = 20, quite above the critical value so that analysis on the 
low-temperature phase becomes easier. 

The first observation is made for the helicity modulus. As shown in Fig. 1, in-plane hclicity moduli set up 
at T ~ 0.96 J/ks- While the difference between the transition points of the two anisotropy parameters is very 
small to be seen clearly by the present system size, we differentiate them as T m and Tkt for 7 = 8 and 7 = 20, 
respectively, for the reasons elucidated later. The finite helicity modulus T y along the magnetic field indicates 
the phase coherence at low temperatures, corresponding to the breaking of U(l) gauge symmetry and therefore 
appearance of superconductivity. A finite value of T x is a reflection of the intrinsic layer pinning (that along the 
c axis remains vanishing down to zero temperature). The collapse of the helicity moduli in x and y directions for 
7 = 20 is an intrinsic property of the state the system presumes beyond the multicritical point. It is considered 
to be responsible to the experimentally observed orientation independence of I-V characteristics, noting that the 
conductance can be evaluated by helicity moduli. 

The second observation is on the density correlations of Josephson vortices. Structure factors S[k x , y = 0, k c ] at 
the transition points are displayed in Fig. 2. As seen in Fig. 2(a), six Bragg peaks appear for 7 = 8 (see also Fig. 
3), indicating that a 3D LR crystalline order is established upon the first-order freezing transition (T m refers to the 
melting point). In a sharp contrast, the Bragg spots at [±2fn/d, ±n/d] for 7 = 20 shown in Fig. 2(b) are diffusive 
and stripe like. As in Fig. 3, the k c profile of the Bragg spots for 7 = 20 is fitted well by a Lorentzian function 
i/[?tT 2 + (k c — Tr/d) 2 ] with the correlation length £ c ~ 0A7d in the c direction. Josephson vortices are therefore 
decoupled into nearly independent layers at the transition point. The correlation length grows to £ c ~ 1.5d at 
T = 0.7 J /ks- The SR crystalline order in the c axis makes the interlayer shear modulus vanishing, which was first 
discussed by Blatter et al. [6]. 

The k x profiles of Bragg spots at [k x , k c ] = [±2fir/d, ±ir/d] for 7 = 20 are plotted in Fig. 4 at several typical 
temperatures. Singularity appears clearly in the k x profile below Tkt = 0.96J//cs, consistent with the onset of the 
hclicity moduli in Fig. 1. By fitting the profiles of the Bragg spots to the function / ~ a\l — fc x d/2/7r| J)_2 + b\l — 
k x d/2fn\ + c, denoted in Fig. 4 by the solid curves, we estimate the exponent as 77 ~ 2.40 ± 0.03, 2.21 ± 0.01 and 
1.87 ± 0.01, 1.41 ± 0.02 for T = 0.96, 0.95 and 0.8, 0.7 J/k B , with the error bars from the least-squares fittings; for 
T = 0.92 J/ks a logarithmic function fits best with the data, corresponding to 77 = 2; data for T = 0.7 J/ks are not 
included for the sake of clarity. The power exponent 77 depends on temperature, implying that elastic constants 
are not simply entropy dominated as in polymer systems. The height of Bragg spots increases very slowly with 
decreasing temperature (see Fig. 1 of Ref.[18] for comparison). The Bragg spots become very sharp to be fitted 
by a power-law function for T < T x ~ 0.65 J/ ks; there is no difference between the shape of structure factors for 
7 = 20 at these low temperatures and that for 7 — 8 in Fig. 2(a). These observations indicate a 3D LR crystalline 
order realized even in highly anisotropic systems provided temperature is low enough. 

The singularities in the structure factors for 7 = 20 at T x < T < Tkt indicate unambiguously power-law 
density correlations, characterized by the exponent 77, along the x direction in real space. Combining with the 
similar observations along the direction of the magnetic field, with smaller exponents, a 2D QLR crystalline order 
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is concluded at these intermediate temperatures. Since all 2D QLROs associated with continuous degrees of freedom 
known to date are governed by the KT fixed point, we identify the intermediate temperature region as a novel KT 
phase. The melting of this 2D Josephson vortex lattice at Tkt is therefore considered as a KT transition, at which 
the tilt modulus diminishes to zero. Since the crystalline order is SR in c direction in the KT phase, so do the 
gauge invariant phase correlations. The finite hclicity moduli for 7 = 20 below T K t correspond to QLROs of phase 
variables. 

The third observation is on the trajectories of the Josephson flux lines. This is important since a decoupling 
between ordered layers is argued based on RG analysis to be impossible, providing Josephson flux lines are confined 
completely by the Cu02 layers [15]. As depicted in Fig. 5, hops of Josephson flux- lines segments across Cu02 layers, 
by creating pancake vortices, into neighboring reservoir layers are observed for T > 0.6J/fcs, and the percentage 
of hopping Josephson flux lines increases quite sharply with temperature. These observations move away the 
hurdle to the decoupling and the 2D phase, and thus reconcile our simulation results with the RG analysis. Hops 
of Josephson flux-line segments disturb intralayer correlations, weaken the interlayer correlations, as clearly seen 
in Fig. 4, and result in the suppression of 3D crystalline order into decoupled 2D quasi lattices for weak bare 
Josephson couplings. For T < T x , hops of Josephson flux lines are suppressed, where the argument by Mikhccv 
and Kolomeisky should apply. As a matter of fact, this temperature regime coincides roughly with that in which 
the 3D LR crystalline order becomes stable for 7 = 20. In Fig. 5, we also display the temperature dependence of 
populations of thermally excited vortex loops, and of collisions between Josephson flux lines residing in the same 
block layers. These results indicate clearly the importance of thermal fluctuations in the present system, and the 
former explains why the exponent r\ depends on temperature. 

The evolution from the 2D QLRO to 3D LRO at T x is also believed to be a thermodynamic phase transition. 
Resorting to the effective Landau theory formulated first by Balents and Nelson [11], which works better for the 
present QLRO in the KT phase, it is argued that this phase transition is probably second order and in the 3D 
XY universality class. This observation can be taken consistent with the absence of noticeable anomaly in the 
specific heat at T x [18, 19], since the critical exponent a is slightly negative in the universality class and thus the 
associated cusp could be very small. The interlayer shear modulus is suppressed continuously as Cqq <~ (d/f£ c ) 2 
with an exponent 2v ~ 4/3 as T x is approached from below. 

Possible size effects on our simulation results are addressed as follows. The successful observation on the single, 
first-order melting transition for 7 = 8 indicates that the system size is sufficient for deriving right physics below 
the multicritical point. Since increasing the anisotropy parameter only reduces the coupling in the c axis, the 
system size unlikely becomes insufficient above the multicritical point. It is easy to see that the 2D QLRO phase 
[T>< , Tkt] does not shrink to zero in the thermodynamic limit, since T KT > T£| ro > T x , with T^ ro ~ 0.89 J/k B , is 
observed under periodic boundary conditions. 

Based on the analyses presented so long, we map out in Fig. 6 the B — T phase diagram for interlayer Josephson 
vortices, noting that the same physics should occur when the magnetic field is tuned while the anisotropy parameter 
is fixed. The first-order melting line for low magnetic fields branches into two phase boundaries at the multicritical 
point B mc = 4>o/2\/3jd 2 , containing an intermediate KT phase. 

The present results are discussed in consistency with previous researches in literature: The KT phase and hops 
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of Josephson flux-line segments observed in the present simulations provide a clear support to the scenario by 
Chakravarty et al. and Blatter et al. [5, 6] formulated in order to explain the peculiar orientation independent I-V 
characteristics and the power-law, non-Ohmic dissipations [3, 4]. The recent experimental observation by Schilling 
et al. [20] on first-order melting of Josephson flux-line lattice up to 10T for YBa 2 Cu 3 07_,5 is consistent with 
B nlc ~ 50T for 7 = 8. The steep normal to superconductivity phase boundary at high magnetic fields observed by 
Lundqvist et al. [21] is able to be explained by the lower bound T^ e on Tkt(B). Ooi and Hirata [22] detected a 
phase boundary on which the 3D triangular lattice softens, very similar to the one in our proposed phase diagram. 

Simulations are performed on the Numerical Materials Simulator (SX-5) of NIMS. This study is partially 
supported by MEXT, Japan, under the Priority Grant No. 14038240. 
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Figure captions: 

Fig. 1: Temperature dependence of in-plane hclicity moduli for 7 = 8 and 7 = 20. The universal drop of helicity 
modulus for a pure 2D XY model at the bare KT transition is expected at T ~ 0.56 J and T^ c ~ 0.89 J/ks- 
Fig. 2: Structure factors at the transition point T^kt = 0.96J/fcs for 7 = 8 (a) and 7 = 20 (b). 
Fig. 3: k c profiles of the Bragg peaks at [k Xl k c ] = [±2fir/d,±ir/d] in Fig. 2. The solid curve for 7 = 20 is the 
result of the least-squares fitting to the Lorentzian function as described in text. 

Fig. 4: k x profiles of the Bragg spots at [k x ,k c ] = [±2fn/d, ±ir/d] for 7 = 20 at several typical temperatures. The 
solid curves are results of the least-squares fittings to the power-law function as described in text. 
Fig. 5: Temperature dependence of ratios of Josephson flux lines which contain segments hopping into neighboring 
block layers, of those which hop and/or collide with neighbors in the same block layers, and populations of thermally 
excited, closed loops of Josephon vortices (normalized by 40 x 240) and those containing also pancake vortices 
(normalized by 240). 

Fig. 6: B — T phase diagram for interlayer Josephson vortices with a multicritical point. The phase boundaries 
T m (B), Tkt(B) and T X (B) are associated with first-order, KT and 3D XY phase transitions as discussed in text. 
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